A texture tensor to quantify deformations: the example of two-dimensional flowing 

foams 
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Abstract: In a continuum description of materials, 
the stress tensor field 7r quantifies the internal forces the 
neighbouring regions exert on a region of the material. 
The classical theory of elastic solids assumes that ef de- 
termines the strain, while hydrodynamics assumes that 
a determines the strain rate. To extend both successful 
theories to more general materials, which display both 
elastic and fluid properties, we recently introduced a de- 
scriptor generalizing the classical strain to include plastic 
deformations: the "statistical strain", based on averages 
on microscopic details ( "A texture tensor to quantify 
deformations" M.Au., Y.J., J.A.G., F.G, companion pa- 
per, Granular Matter, same issue). Here, we apply such 
a statistical analysis to a two-dimensional foam steadily 
flowing through a constriction, a problem beyond reach 
of both theories, and prove that the foam has the elastic 
properties of a (linear and isotropic) continuous medium. 



A "plastic" deformation means that microscopic rear- 
rangements take place in the material, so that the mi- 
croscopic pattern does not return to its initial condition 
even after the applied force has ended. An example of 
a highly heterogeneous one in a viscoelastic material is 
a two-dimensional foam steadily flowing through a con- 
striction (Fig. 1). This apparently simple example is 
utterly intractable from the perspective of both elastic- 
ity theory [1] and Navier-Stokes [2] treatments. In this 
paper, we present a new approach to analyse complex 
flows of disordered materials. 

We prepare the foam by blowing air into the bottom of 
a column of soap solution. The solution is 10% commer- 
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cial dishwashing detergent ("Ivory" brand) and 5% glyc- 
erol. Its surface tension, measured by the Noiiy method, 
is 7 = 28.5 ±0.1 mN.ni -1 . Filtered air blows at a steady 
flow rate of 0.08 cm 3 .s _1 through an 18-gauge (.084 cm 
ID) stainless steel needle with a 90° bevel. Bubbles are 
homogeneous in size (<5% dispersity), and float to the 
top of the column in random positions. The foam is dry, 
with a relative fraction of fluid <3%. Since a bubble edge 
(a film of soap solution) has two interfaces with air, its 
line tension A is A = 2~fh = 28.5 //N. 

Bubbles enter a horizontal channel, a Hclc-Shaw cell 
made of two parallel Plcxiglas plates h — 0.5 mm apart, 1 
m long, and 10 cm wide (Fig. 1). The average bubble ve- 
locity is 0.1 cms -1 . The generation of new bubbles forces 
the mass of bubbles down the channel, in plug flow (free 
slip boundary conditions) if the channel was uniform. A 
5-mm wide constriction near the end of the channel forces 
the flow to be heterogeneous. Beyond the constriction the 
bubbles reexpand into the full width channel for a short 
distance and then fall from the open channel end (no 
stress boundary conditions) into a collection vessel. The 
volume variation during the transit time (50 s) across 
the field of view, due to pressure differences or diffusion 
of gas from one bubble to its neighbours, are below our 
pixel induced detection limit of 2%. 

FIG. 1. Two-dimensional foam flowing through a constric- 
tion. The 10 cm (551 pixels) wide field of view shows only 
the end of the 1 m long horizontal channel. 

All measurements we present use 60x60 pixels sliding 
boxes ("representative volumes"), meaning that we do 
not evaluate them within 30 pixels of the channel walls 
; and average over 2800 successive images of a 30 Hz 
movie. To measure the "Eulerian" velocity field, we track 
the center of mass of each bubble between two successive 
images and add the velocities of all bubble in the same 
volume. The velocity field is smooth and regular (Fig. 
2), qualitatively indicating that the foam behaves as a 
continuous medium. 

FIG. 2. Velocity field in the foam (in arbitrary units, the 
same scale for each arrow). 

To obtain a more quantitative characterization, we 
measure the stress in the foam. Stress has dissipativc 
and elastic components; the pressure inside the bubbles 
and the network of bubble edges contribute to the latter. 
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Since the pressure stress is isotropic, it does not con- 
tribute to the elastic normal stress difference a xx — a yy 
(as measured for granular materials [3]) or shear stress 
W xy , entirely due to the network. Dimensionally, a is of 
the order of the line tension A (which, in a 2D foam, is 
the same for all edges: here 28.5 /zN) divided by a typical 
bubble size. 

We measure locally the network stress in each "repre- 
sentative volume element" that is a square box, centered 
around the point of measurement, of a mcsoscopic size: 
larger than a bubble, but much smaller than the channel 
width. We proceed as follows. We identify the bub- 
ble edges which cross the boundaries of the volume [4]. 
We determine the tension f = Ae of each edge, where 
e is the unit vector tangent to the edge. We determine 
the average force / on a boundary element dS by vec- 
torially adding these tensions and obtain a defined by: 
fi = aijdSj [1,5]; equivalently, we use an average over all 
links to improve the statistics [6,7]. 

Clearly, the stress field is strongly heterogeneous (Fig. 
3). The upstream influence of the constriction becomes 
visible as the lobe where a xx — W yy changes sign. 

FIG. 3. Snapshot of the foam (Fig. 1) with regions colored 
according to the sign of the experimental value of the normal 
elastic stress difference a xx — a yy : blue x, negative; black o, 
zero within error; red +, positive; green o, values we omit in 
Fig. 5. 

To characterize the deformation of the bubbles, we con- 
sider a volume V(R) around the position R. In this vol- 
ume, we list all pairs (r, r 1 ) of positions of neighboring 
vertices connected by one bubble edge. From the vector 
f = r'-f, we construct the tensor I® I = This ten- 

sor averaged over V(R) defines the local "texture tensor" 
W(R): 

see details in the companion paper [7]. This symmet- 
ric tensor has two strictly positive eigenvalues, both of 
order {t 2 ), the largest being in the direction in which 
bubbles elongate. It reflects at large scales the relevant 
features of the actual microstructure of_the material. For 
instance, Fig. 4 shows an example of M(R) measured in 
the flowing foam experiment of Fig. 3. This texture ten- 
sor M(R) thus quantifies the qualitative impression of 
compression or elongation we obtain by looking at the 
bubbles in each region of Fig. 1. The smooth and con- 
tinuous variations of the statistical strain tensor field in 
space validate treating the foam as a continuous medium. 



FIG. 4. Experimental measurement of M(R), superim- 
posed on a snapshot of the foam (Fig. 1). Since M is a 
real symmetric tensor with strictly positive eigenvalues, we 
display it as an ellipse, the length of each axis (in arbitrary 
units, the same scale for each ellipse) being proportional to 
an eigenvalue. We omit data on boxes touching the channel 
wall. 

We have defined the "statistical strain tensor" [7] : 
V{R) = i (\ogM{R) - logM ) . 

Here the reference value M is chosen in the undeformed, 
isotropic foam far upstream of the constriction (a choice 

that plays no role in elastic properties): M = (£l)I/2, 
where / is the identity tensor, here in two dimension. 
By definition, we calculate the tensor logM by rotating 
M to a basis where it is diagonal, take the logarithm of 
its eigenvalues, then rotate it back into the initial basis 
(this definition is more convenient than the equivalent 
one using a development as an infinite scries). Hence 

log M has the same axes as M, and is real and symmetric. 

The "statistical strain" tensor U reduces to the usual 
definition of strain u in the validity limits of classical 
elasticity [7]. It is purely geometric, and does not explic- 
itly depend on stresses and forces. It applies to a whole 
general class of materials with both elasticity and plastic 
rearrangements, whether 2D or 3D, whenever we can ex- 
perimentally measure the relevant information: whether 
detailed (list of microscopic positions and links) or meso- 

scopic averages (quantities related to M). It is a state 
variable, constant in a steady flow, and is not necessar- 
ily homogeneous through the whole sample, allowing a 
thermodynamic description of non-equilibrium complex 
fluids [8]. 

Do the elastic stress and the statistical strain relate 
deterministically? In Fig. 5 we plot the normal stress 

difference a xx — a yy versus U xx — U yy ; in this represen- 
tation, the isotropic contributions (of Mq in U, and of 
the pressure stress in a) play no role, hence we refer only 
to the foam's current state. Each data point is a mea- 
surement which derives from averages at one position of 
the foam: data sources are the same as in Figs. 3 and 
4. That is, in principle Fig. 5 measures a constitutive 
relation from a single image; in practice, we average 2800 
similar images to improve the precision. 

Since M, and hence U, is completely independent of 
<7 [7] , the high correlation between U and a apparent in 
Fig. 5 reflects the physical constitutive relation required 
to treat the foam as a continuous medium, in which de- 
tails of the microstructure appear only through meso- 
scopic averages. Since different applied strains u apP i can 
correspond to the same a, such relation does not appear 
in classical a vs. u app i plots [6]. 
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FIG. 5. Plot of U versus a. The color code is the same 
as in Fig. 3, showing the sign of the normal stress difference 
~<Jxx — ~<Jyy'- blue, negative; black, zero within error bars; red, 
positive. We omit data from the boxes touching the channel 
wall. The black line is a linear fit. On the same graph, we plot 
a xy (violet) versus U xy , shifted by —1 mN.m -1 for clarity. 

Moreover, the relation between U and cf is linear over 
the whole range covered by this experiment. The slope 
of Fig. 5 measures 2fi, where fi is the shear modulus of 
the foam, much beyond the validity regime of classical 
elasticity. We find, for the xx — yy component: 

pt xx _ vy = 3.00 ± 0.005 mN.rrT 1 . 

The foam is nearly isotropic, since we find almost the 
same value for the xy component: 

fj, xy = 3.06 ±0.01 mN.rrT 1 . 

Note that we find similar values, 3.00 ± 0.01 and 3.01 ± 
0.01 respectively, for a foam flowing relatively to a fixed 
round obstacle (data not shown). A regular "honey- 
comb" hexagonal lattice has an isotropic shear modu- 
lus Hh — X/(£y/3), hence with the same tension and av- 
erage bubble area (A) it would have a shear modulus 
Hh = 0.465 A/VW = 2.52 mN.nT 1 [9-11]. Our value 
/i = 3.0 is higher: /J,//J,h = 1-2. This difference is likely 
due to the side length disorder (i 2 ) / {if = 1.22 of the 
foam. How fi/fih depends on the microstructure and its 
disorder is an open problem [12]. 

We have performed the following tests of the validity of 
our method. We have checked the agreement between our 
measurement of statistical strain and the classical strain, 
first analytically on a honeycomb lattice, then numeri- 
cally on a disordered simulated foam kept in the elastic 
regime. We have checked that if we increase the size of 
the volume element, we decrease the fluctuations, and 
thus the upper and lower values of a and U in Fig. 5, 
but the slope and hence /i are unchanged. Finally, the 
value of the shear stress tensor field or shear modulus, 
resulting from our image analysis, predicts the force ex- 
erted by the flowing foam on an obstacle. We have tested 
it against an independent, macroscopic measurement of 
the force, and found that they agree (S. Courty et al, in 
preparation) . 

In summary, we have measured the stress, the texture 
tensor and the statistical strain for a 2D flowing foam. 
We have shown that the 2D foam behaves like a linear 
and isotropic continuous material. 
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This figure "Fig3artExp.gif" is available in "gif" format from: 
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